AD-AU'7  9*+2 

UNCLASSIFIED 


OKLAHOMA  UNI V  NORMAN  SCHOOL  OF  AEROSPACE  MECHANICAL  —ETC  F/S  1 
MECHANICS  OF  BIMODULAR  COMPOSITE  STRUCTURES. (U) 

J  N  REDDY  N00018-78-C-0647 


I 


Department  of  the  Navy 
OFFICE  OF  NAVAL  RESEARCH 
Mechanics  Division 
Arlington,  Virginia  22217 

Contract  N00014-78-C-0647 
Project  NR  064-609 
Technical  Report  No.  28 


Report  OU-AMNE-82-4  &  VPI-E-82 .20 


MECHANICS  OF  BIMODULAR  COMPOSITE  STRUCTURES 


C.W.  Bert 

School  of  Aerospace,  Mechanical  and  Nuclear  Engineering 
University  of  Oklahoma 
Norman,  Oklahoma,  USA  73019 


J.N.  Reddy 

Department  of  Engineering  Science  and  Mechanics 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia,  USA  24061 


July  1982 


University  of  Oklahoma 
Norman,  Oklahoma  73019 


DTIC 

\>"  AUG  0  61982 


Approved  for  public  release;  distribution  unlimited 


82  08  06  nn 


Advances  in  Mechanics  of  Composite  Materials 

(Proceedings  of  che  I  IT  AM  Symposium,  Virginia 
Polycech.  Insc.,  Blacksburg,  VA,  Aug.  1982) 


MECHANICS  OF  BIMODULAR  COMPOSITE  STRUCTURES 
C.W.  3erc*  and  J.N.  Reddv** 


V; 

O  \i  I 


I'  i  r;  (.  r 1 
Av-  i  1  ■ 


Dir.t 


a 


This  paper 
fiber-rein 
properties 
cord-rubbe 
Specific  c 
materials 
inaced  and 
tr ansienc 
thickness- 
methods  us 


‘School  of  Aerospace,  Mechanical  and  Nuclear  Eng., 

Untv.  of  Oklahoma,  Xcrr.cn ,  Oklahoma,  USA  73019 
**Dcpt.  of  Eng.  Sci.  and  Mechanics,  Virginia  Polytech.  Inst, 
and  State  Univ. ,  Blacksburg,  Virginia,  USA  -1061 


ABSTRACT 

is  a  survey  of  the  mechanics  of  bean  and  plate  structures  laminated  of 
forced  composite  materials  having  different  elastic  and  theraoelastic 
in  tension  and  compression.  Examples  of  such  materials  include  tire 
r,  wire-reinforced  solid  propellants,  and  soft  biological  materials, 
opics  covered  include:  mathematical  models  of  fiber-reinforced  bimodulcr 
and  their  exoerimental  verification;  static  and  dynamic  analysis  of  lam- 
sandwich  beams ;  plane  elasticity;  analysis  of  deflection  and  free  and 
vibration  of  laminated  plates  and  shells.  In  all  of  these  analyses, 
shear  deformacion  and  rotatory  inertia  are  included.  The  solution 
ed  include  closed-form,  transfer-matrix,  and  finite-element  techniques. 
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INTRODUCTION 

The  knowledge  that  certain  materials  have  different  stress-strain  behavior  when 
they  are  loaded  in  tension  in  contrast  to  compression  goes  back  more  than  a  cen¬ 
tury  (Saint -Venant ,  1864).  However,  apparently  Timoshenko  (1941)  was  che  first  to 
consider  che  bilinear  approximation  for  che  uniaxial  stress  state,  with  one  slope 
(modulus)  in  censior.  and  a  distinctly  different  one  in  compression.  Modem  inter¬ 
est  in  the  topic  stems  from  the  general  isotropic  formulation  of  Ambartsumyan 
(1965),  who  is  credited  with  originating  the  "bimodulus"  or  "bimodular"  terminology. 

Experimental  evidence  of  drastically  different  stress-strain  behavior  in  tension 
and  compression  has  been  found  in  numerous  materials,  including  cast  iron  (Gilbert, 
1961),  tire  cord-rubber  (Clark,  1963),  polyamides  (Zemlyakov,  1965),  polycrystal¬ 
line  graphites  (Seldin,  1966),  wire-reinforced  solid  prooeilants  (Herrmann  et  al., 
1967),  concrete  (Seefried  et  al.,  1967),  epoxies  (Novak  and  3ert ,  1963),  carbon- 
carbon  composites  (Adsic  et  al.,  1972),  cortical  bone  (Simkin  and  Robin,  1973), 


2 


1974),  compacted  clay  soils  (Ajaz  and  Parry,  1975),  sintered  metal  (Ducheyne  et 
al. ,  1978),  soft  biological  tissues  (Pearsall  and  Roberts,  1978),  and  paperboard 
(Carlsson  et  al.,  1980);  see  Table  1.  It  should  be  remarked  that  the  no-tension 
material  model  often  used  in  rock  mechanics  (Zienkiewicz  cl  al.,  1968)  is  an 
extreme  limiting  case  of  bimodular  material.  For  an  extensive  survey  of  bimcdular 
material  behavior  as  well  as  the  strength-differential  phenomenon  (different 
strengths  in  tension  and  in  compression),  the  reader  is  referred  to  the  work  of 
Kamiya  (1976a) . 

TABLE  1  Some  Examples  of  Differpr.t  Elastic  Behavior  in 

_ Tension  and  Compression _ 


Investigator 

Material 

Young's  Moduli  Ratio 
EC/Et 

Kotlarskii  6  Karbasova  (1968) 

Fabric/ rubber 

0.38 

Ducheyne  et  al.  (1978) 

Sintered,  porous 
stainless  steel 

0.1 

Zolotukhina  &  Lepetov  (1968) 

Various 

fabric/rubber 

0.07  to  0.50 

Patel  et  al.  (1976) 

Polyester 

cord/rubber 

0.017 

Aramid 

cord/rubber 

0.0034 

Posfalvi  (1977) 

Rayon 

cord/rubber 

0.0036 

Pearsall  Roberts  (1978) 

Myometrium 

(uterine  muscle) 

0.2 

Here,  initial  emphasis  is  placed  on  the  mechanics  of  fiber-reinforced  composites 
with  bimodular  action  on  two  different  levels:  micromechanics  to  explain  the 
physical  phenomena  involved  and  anisotropic  elastic  continuum  theory  for  use  in 
structural  mechanics.  This  is  followed  by  applications  to  laminated  and  sandwich 
beams,  plane  elasticity,  plate  bending,  and  shells. 


MI CROMECltAN  1 CS  OF  COMPOSITES  WITH  CURVED  FIBERS  AND  SOFT  MATRICES 

To  illustrate  the  severity  in  bimodular  material  action,  we  present  stress-strain 
data  in  Fig.  1  for  unidirectional  polyester-cord-reinforced  rubber,  taken  from  the 
work  of  Bert  and  Kumar  (1981).  The  basic  phenomena  taking  place  are  fiber 
straightening  due  to  tie-bar  action  in  tension  and  fiber  curving  due  to  column 
action  in  compression.  Both  of  these  phenomena  are  reacted  or  restrained  by 
interaction  of  the  fibers  with  the  matrix,  which  acts  like  an  elastic  foundation. 
The  first  model  based  on  these  phenomena  is  due  to  Herrmann  et  al.  (1967).  A  more 
recent  micromechanics  model  is  due  to  Cominou  (1976).  Numerous  models  of  this  type 
were  mentioned  in  a  recent  paper  (Bert,  1979),  in  which  the  effects  of  fiber  ten¬ 
sion  or  compression  on  thermal  expansion  coefficient  were  also  predicted.  The 
latter  are  in  qualitative  agreement  with  the  experimental  thermal-expansion 
results  of  Makarov  and  Nikolaev  (1971). 

Another  analytical  approach  to  the  micromechanics  of  composites  with  curved  fibers 
is  the  mean-fiber-angle  concept  originated  by  Tarnopol ' skil  et  al.  (1967).  They 
assumed  that  the  fibers  are  initially  sinusoidally  curved  with  an  initial  deflec¬ 
tion  given  by 

w  -  VI  sin  8  x 
o  o 


(1) 


3 


Fig.  1.  Stress-strain  curve  for  polyester-rubber  in  the 
cord  direction  (1  psi  =  6.895  x  103  pascals). 

Then  the  inclination  at  an  arbitrary  axial  position  x  is 

6  -  arctan  (6  W  cos  3  x)  (2) 

o  o  o 

Using  equation  (2)  in  the  veil-known  transformation  equations  for  Young's  modulus 
and  integrating  over  a  half-wave  length  2,  they  obtained  the  following  result  for 
the  f iber-direccion  initial  Young's  modulus  of  a  wavy-fiber-reinf orced  composite: 

a/Ex)  =  (l/2Ej)  (2  +  f02)(l  +  f0V3/2 

+  (G12-X  -  2v12El“l>(f02/2)(l  +  f0V3/2 

+  (1/2E_)  [  2  -  (2  +  3f  z)(l  +  f  2)'3/2]  (3) 

2  o  o 

Here  and  E2  are  the  major  (fiber  direction)  and  minor  Young's  moduli  of  a  uni¬ 
directional  composite  with  straight  fibers,  E^  is  the  major  Young's  modulus  of  the 
vavy-fiber  composite,  f  =  B  W  ,  and  G,,  and  v.,  are  the  respective  major  shear 

O  O  O  x.c  14 

modulus  and  major  Poisson's  ratio  of  the  straight-fiber  composite. 

For  f  <<  1,  the  following  approximate  expression  was  suggested 

E1/E1  :  [1  +  (fo2/2)(E1/C12))"1  (4) 

The  disadvantage  of  using  either  equation  (3)  or  (4)  is  that  it  requires  prior 
knowledge  of  a  straight-fiber  composite  having  the  same  fiber  volume  fraction. 

For  the  glass-fiber  reinforced  plastic  investigated  experimentally  by  Tarnopol ' skii 
et  al.  (1967),  it  was  found  empirically  that  an  initial  composite  prestress  of 
approximately  107;  of  the  ultimate  tensile  strength  was  necessary  to  assure  that 
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the  fibers  are  sufficiently  straight.  This  work  can  be  used  to  predict  only  the 
Initial  modulus  (or  Poisson's  ratio).  However,  in  certain  composites  wich  a  highly 
flexible  matrix  (such  as  rubber),  the  composite  modulus  may  change  quite  signifi¬ 
cantly  with  loading  (especially  in  compression,  as  discussed  earlier  in  chic  sec¬ 
tion).  Thus,  the  analysis  of  Tabaddor  and  Chen  (1971)  is  quite  Important  for  such 
composites,  since  it  permits  computation  of  the  complete  stress-strain  curve  in 
both  tension  and  compression  for  the  composite  in  the  fiber  direction.  In  their 
work,  they  replaced  f  in  equation  (3)  with  f  (  =  3W)  and  developed  the  following 

relationship  between  f  and  fQ,  where  W  is  defined  in 

w  •  W  sin  Sx  (5) 

f  -  BW  -  3  W (l  +  c  )/(l  +  t  )  (6) 

o  o  i  1 


Later  Tamopol'skii  et  al.  (1973)  investigated  the  effect  of  fiber  waviness  in 
composites  with  anisotropic  fibers  (specifically  carbon  fibers).  They  found  that 
for  such  composites,  it  is  necessary  to  use  the  full  equation  (3)  rather  than  the 
Simplified  equation  (9).  Van  Dreumel  and  Kamp  (1977)  investigated  experimentally 
the  effects  of  fibei  waviness  on  both  tensile  strength  and  major  Poisson's  ratio 
of  carbon-fiber  reinforced  plastic. 


In  order  to  completely  characterize  a  thin  sheet  material  for  use  in  stress  analy¬ 
sis  and  structural  design,  one  needs  a  complete  set  of  elastic  stress-strain 
relations,  for  fiber-direction  tension  and  compression  (see  the  next  section),  not 
just  the  Young's  moduli.  A  rational  analytical  basis  for  such  a  complete  charac¬ 
terization  was  developed  by  Bert  (1979).  His  expressions  for  the  curved-fiber 
composite  compliances  S  in  terms  of  the  straight-fiber  unidirectional  compliances 
are  as  follows:  J 


fhJ 

S22 

X1 

I2_ 

'  r 

!“ 

►  = 

S12 

X2 

-1! 

7 

S22 

S11 

*3 

Z2 

.  S&6 

*4 

-X4_ 

v  *  > 

Here  the  invariants  are 


h  S  2S12 


2S22  +  S66 


I3  ,  -  2SU  +  2S12  +  Sfe6 


I2  5  S11  "  2S12  +  S22  ~  S66 
*4  H  4S11  -  8S12  +  4S22  -  4S66 


(8) 


The  average  values  of  the  functions  of  m  (  =  cos  8)  and  n  (=  sin  0)  appearing  in 
equations  (7)  are 


_  L/4  L/4 

»4  =  [2  [  m4  dx) /[ 2  f  dx]  -  [1  +  ( f 2 /2) ) (1  +  f2)_3/2 
J-L/4  J-L/4 

T  =  (1  +  f2)-1/2  ,  m3n  -  mn3  -  0 

D 


(9) 


3  3 

It  is  noted  that  since  m  n  ■*  mn  »  0,  the  wavy-fiber  composite  remains  orthotropic. 
Use  of  equations  (7)-(8)  in  an  iterative  manner  permits  development  of  complete 
nonlinear  strain  versus  stress  relationships  for  the  composite. 
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The  experiments  of  Patel  et  al.  (1976)  showed  that  the  difference  between  tension 
and  compression  properties  were  much  more  pronounced  for  rubber  reinforced  with 
aramid  cord  than  when  reinforced  with  polyester  cord.  This  is  probably  related 
to  the  very  low  compressive  strength  of  the  aramid  fibers  themselves,  which 
Greenwood  and  Rose  (1974)  suggested  is  duo  to  separation  of  the  xicroflbrils  com¬ 
prising  the  fibers.  Moncunill  de  Eerran  and  Harris  (1970)  found  experimentally 
that  the  buckling  mode  of  steel  wire/polyester  was  helical  rather  than  planar  as 
assumed  in  all  known  analyses.  It  has  been  suggested  that  the  more  pronounced 
difference  between  tensile  and  compressive  behavior  exhibited  in  the  tire-cord/ 
rubber  composites  of  Patel  et  al.  (1976)  ard  Posfalvi  (1977)  as  compared  to  that 
of  straight  wire-reinforced  rubber  may  be  due  to  the  inherent  twisting  action  pre¬ 
sent  in  tire  cord,  which  has  a  helical  geometry. 


MACROMECHANICS  OF  Bt MODULAR  ANISOTROPIC  CONSTITUTIVE  RELATIONS 

The  first  work  on  the  mechanics  of  bimodular  orthotropic  media  is  due  to 
Ambartsumyan  (1965,  1969).  Other  early  models  are  due  to  Isabekyan  and 
Khachatryan  (1969)  and  Shapiro  (1966).  More  recent  work  in  this  area  intended 
especially  for  transversely  isotropic  ATJ-S  graphite,  is  due  to  Jones  and  Nelson 
(1976),  Jones  (1977a, b),  and  Jones  and  Morgan  (1977). 

For  fiber-reinforced  composites  with  very  soft  matrices  (e.g,,  .ord-rubber) ,  Bert 
(1977)  introduced  Che  .iber-govemed,  symmetric  compliance  model  and  applied  it 
very  successfully  to  experimental  data  reported  by  Patel  et  al.  (1976).  In  this 
model,  chere  are  two  sets  of  elastic  coefficients:  one  when  the  fiber-direction 
strain  (e^)  is  tensile  and  another  when  is  compressive.  Within  each  set 

(tensile  and  compressive),  the  compliance  matrix  is  symmetric,  so  that  the  concept 
of  strain  energy  is  retained.  Thus,  within  this  context,  the  well-known  reciprocal 
relation  relating  Young's  moduli  (Ej)  and  Poisson's  ratios  (v^j)  holds: 

vLT/EL  "  VTL/ET  (k«c,c)  (10) 


Here,  set-scripts  L  and  T  refer  to  the  directions  respectively  parallel  with  and 
perpendicular  to  the  fiber  direction,  and  superscripts  c  and  t  refer  to  fiber- 
direction  compression  or  tension  strain.  It  should  be  emphasized  that  E£  and 


although  referred  to  as  compressive  properties  according  to  the  convention  men¬ 
tioned,  are  actually  obtained  under  uniaxial  tensile  loading  perpendicular  to  the 
fibers.  In  particular,  the  four  properties  obtained  from  uniaxial  tensile  tests, 
t  t  c  c 

namely  E^,  v^,  E^,  and  v.^,  are  not  related. 


A  linear  material  model  can  be  characterized  by  either  Its  compliance  matrix  or 
alternatively  by  its  stiffness  matrix.  The  criteria  used  to  evaluate  the  various 
bimodular  material  models  are  as  follows: 

1.  The  model  should  be  related  to  the  governing  physical  mechanisms 
involved,  i.e.,  fiber  tie-bar  effect  in  tension  and  fiber  buckling  instability  in 
compression. 

2.  It  is  highly  desirable  to  have  the  strain  energy  be  positive 
definite  so  that  energy  is  conserved. 

3.  The  compliances  should  be  consistent  with  measured  values  for  the 
conditions  specified,  i.e.,  they  should  depend  upon  the  nature  of  the  mulciaxial 
stress  or  strain  state  as  appropriate. 

4.  The  shear  modulus  for  an  orthotropic  bimodular  material  should  have 
different  values  for  shear  stresses  of  same  magnitude  but  opposite  sign  In  any 
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coordinates  other  than  the  material-symmetry  directions. 

As  pointed  out  by  Voigt  (1928),  criterion  2  implies  that:  (a)  the  compliance  or 
stiffness  matrix  be  symmetric  and  (b)  certain  limits  exist  on  the  compliances  or 
stiffnesses  so  that  tne  compliance  or  stiffness  matrix  is  positive  definite. 
Symmetry  of  the  compliance  or  stiffness  matrix  is  necessary  in  ord;r  for  a  material 
to  be  mechanically  stable,  as  shown  by  Srun  (1976).  Furthermore,  compliance  sym¬ 
metry  is  highly  desirable  in  chat  most  structural  analysis  algorithms  are  based  on 
this  assumption,  i.e.,  chey  have  no  provision  for  unsymraetric  compliance  or  stiff¬ 
ness  matrices. 

A  detailed  discussion  and  critique  of  the  various  mathematical  models  for  bi- 
modular  anisotropic  elastic  materials  was  presented  by  Bert  (1978).  Table  2  is 
a  bt.lef  summary  of  the  results  of  this  evaluation  on  the  basis  of  the  four  cri¬ 
teria  mentioned  previously. 

TABLE  2.  Evaluation  of  Various  Macroscopic  Bimodular-Material  Models 


Model 

Model 

Name 

Name 

Cri  terion 

1 

2 

3 

A 

1 

Ambartsumyan  (1965,  1969) 

No 

Yes 

11 

Isabekyan-Khachatryan  (1969)  restricted  compliance 

No 

No 

III 

Shapiro  (1966)  first  invariant 

No 

No 

IV 

Jones  (1977a)  weighted  compliance 

No 

Yes 

Yes 

V 

Bert  (1977)  symmetric  compliance 

Yes 

Yes 

Yes 

Yes 

MECHANICS  OF  BIMODULAR  COMPOSITE-MATERIAL,  LAMINATED,  AND  SANDWICH 
BEAMS 

A  fairly  extensive  survey  of  Che  literature  of  bimodular  beams  was  included  in  the 
recent  work  of  Tran  and  Bert  (1982).  In  the  present  work,  however,  attention  is 
restricted  to  composite  beams.  Many  of  these  were  related  to  the  experimental 
flexural  behavior  of  specific  composite  materials.  Vierling  and  Scheele  (1961) 
were  concerned  with  unidirectional  cord-reinforced  rubber  and  Zolotukhina  and 
Lepetov  (1968)  with  fabric-reinforced  rubber.  Bilida  (1969),  Nachlinger  and 
Leininger  (1969),  Jones  (1976),  Topper  et  al.  (1978),  Zweben  (1978),  and  Phan-Thien 
(1981)  were  involved  with  fiber-reinforced  plastics,  Carlsson  et  al.  (1980)  with 
paperboard,  and  Swift  and  Smith  (1979)  with  reinforced  concrete.  It  should  be  men¬ 
tioned  that  Zweben  (1978)  and  Swift  and  Smith  (1979)  assumed  that  the  materials 
were  linearly  elastic  in  tension  and  perfectly  plastic  in  compression.  The  effects 
of  transverse  shear  deformation  were  included  by  Tabaddor  (1976)  and  Tran  and  Bert 
(1982)  for  static  loading  and  by  Bert  and  Tran  (1982)  for  transient  loading. 

Relatively  little  attention  has  been  devoted  to  the  analysis  of  beems  laminated  of 
bimodular  material.  Kotlyarskii  and  Karbasova  (1968)  were  concerned  with  a  multi¬ 
layer  conveyor  belt,  Crawford  (1968)  with  a  boron-polymer  laminate,  and  Jones  and 
Morgan  (1980)  with  laminates  of  composite  materials.  Recently,  Bert  and  Clifford 
Rebello  (1982)  added  the  effect  of  transverse  shear  deformation  and  made  a  detailed 
Investigation  of  the  effects  of  stacking  sequence. 

The  first  investigation  of  bimodular  sandwich  beams  is  due  to  Bert  (1982),  who 
developed  equations  for  the  reduction  of  experimental  data  on  sandwich  beams  of 
bimodular  material  for  determination  of  separate  stress-strain  curves  In  tension 
and  in  compression  for  the  facing  material.  This  technique  has  a  practical 
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experimental  advantage  over  the  use  of  compact -section  beams  (Nadai,  1950)  In  that 
it  Is  not  necessary  to  measure  the  beam  deflection.  This  technique  was  used  to 
obtain  0*  compression  loading  in  the  work  of  Bert  and  Kumar  (1981),  see  Fig.  1 
of  the  present  paper. 

Apparently  the  only  investigation  of  the  dynamic  behavior  of  sandulch  beams  with 
bimodular  facings  is  due  to  Bert  and  Clarence  Rebello  (1982).  They  presented  both 
classical  (fourth-order  bending  differencial  equation)  and  sixch-order-theory 
analyses  of  pin-ended  beams  and  compared  the  results  with  experimental  results 
obtained  for  beams  with  polyurethane-foam  cores  and  aranid  cord-rubber  facings. 


ORTHOTROPIC,  BIMODULAR  PLANE  ELASTICITY 

Apparently  the  first  formulation  of  the  plane  problem  of  elasticity  for  a  bimodu¬ 
lar,  rectilinearly  orthocropic  material  is  due  to  Ambartsumyan  (1969).  However, 
he  did  not  present  any  solutions.  For  a  bimodular,  cyl indrically  orthotropic 
material,  Kamiya  (1977)  formulated  and  solved  the  plane-strain  thermoelascic  prob¬ 
lem  for  the  axisymmetric  case.  Later  Kamiya  (1979)  generalized  this  work  to 
Include  the  effects  of  temperature-dependent  properties. 

Tabaddor  (1981)  formulated  plane  elasticity  problems  for  rectilinearly  orthotropic 
materials  by  the  finite-element  approach  and  solved  the  case  of  a  point -loaded , 
clamped-c lamped  beam  treated  as  a  plane-stress  problem. 


ORTHOTROPIC  AND  LAMINATED  PLATES  UNDERGOING  BENDING 

Apparently  the  first  formulation  and  solution  of  a  place-bending  problem  for  ortho¬ 
cropic  bimodular  material  was  due  to  Bert  and  Kincannon  (1979).  They  used  Bert's 
fiber-governed  symmetric  material  model  (1977)  and  pointed  out  the  importance  of 
the  neutral-surface  position  associated  with  vanishing  total  strain  in  the  fiber 
direction.  The  particular  problem  that  they  solved  was  a  specially  orthotropic, 
thin  plate  having  an  elliptic  planform,  a  uniformly  distributed  pressure  loading, 
and  clamped  edge. 

The  first  solution  of  a  laminated  bimodular  plate  of  finite  dimensions  was  due  to 
Kincannon  et  al.  (1980) .  They  considered  the  two-layer,  cross-ply  version  of  the 
aforementioned  orthotropic  elliptic  plate  problem.  In  a  cross-ply  laminate,  the 
location  of  the  neutral  surfaces  associated  with  each  of  the  two  fiber  directions 
is  important.  This  is  illustrated  in  Fig.  2,  taken  from  Bert  et  al.  (1981a).  In 
this  figure,  the  neutral-surface  position  (znx)  associated  with  the  bottom  layer 
(fibers  in  the  x  direction)  is  positive (below  the  midplane).  Thus,  compressive 
properties  must  be  used  for  the  upper  portion  (0  <_  z  £  z  )  of  this  layer  and 

tensile  properties  for  the  bottom  portion  (znx  <_  z  <_  h/2,  where  h  5  total  thickness 

of  the  laminate).  The  neutral-surface  position  (zny)  associated  with  vanishing 
total  strain  in  the  y  direction  governs  the  top  layer  and  is  shown  negative,  thus 
falling  within  the  top  layer  and  dividing  it  into  two  regions.  In  summary,  in  this 
most  general  case,  bimodular  action  in  this  two-layer  laminate  causes  it  to  act  as 
if  it  were  a  four-layer  laminate.  Of  course,  It  must  be  remembered  that  the  four 
"layers"  have  different  thicknesses. 

Due  to  their  low  thickness  shear  moduli  relative  to  their  Young's  moduli,  plates 
constructed  of  composite  materials  exhibit  significantly  larger  thickness  shear 
deformations  than  do  isotropic  plates  of  the  same  geometry  and  loading.  Thus,  a 
plate  which  may  reasonably  be  considered  "thin"  (say,  with  a  thickness /s ide  ratio 
less  than  0.1)  for  isotropic  materials  may  have  to  be  considered  "thick"  for 


8 


Fig.  2.  Schematic  diagram  of  an  element  of  a  two-layer,  cross-plv 

plate  constructed  of  bimodular  materials  showing  znx  and  znv, 
the  neutral-surface  positions  associated  with  =  0  and 
ty  «  0.  From  Bert  et  al.  (1981a). 

composite  materials.  The  first  analysis  of  composite-material  plates  including 
thickness  shear  deformation  was  conducted  by  Reddy  and  Bert  (1979).  They  used  che 
finite-element  method  to  analyze  both  single-layer,  specially  orthotropic  and  two- 
layer,  cross-ply  plates  of  rectangular  planform.  The  loading  was  sinusoidally  dis¬ 
tributed  pressure  and  the  boundary  conditions  freely  supported  (simply  supported 
flexurally  with  no  in-plane  normal  restraint)  on  all  four  edges.  Closed-form 
solutions  of  this  problem  were  later  presented  for  thin  plates  and  for  plates  with 
thickness  shear  deformation  (Bert  et  al. ,  1980,  1981a).  Agreement  with  the 
finite-element  results  was  excellent,  thus  providing  a  validation  of  the  finite- 
element  solution. 


The  thermal-stress  problem  of  plates  subjected  to  a  change  in  both  midplane  temper¬ 
ature  and  temperature  gradient  through  the  thickness  was  analyzed  by  Reddy  et  al. 
(1980).  They  presented  both  closed-form  and  finite-element  results  for  rectangular 
plates  with  thickness  shear  and  subjected  to  sinusoidal  distributions  of  midplane 
temperature  and  temperature  gradient. 

It  is  important  to  mention  here  that  in  the  case  of  fiber-reinforced  bimodular 
material,  the  stretching  stiffnesses  (A^),  bending-stretching  coupling  stiffnesses 

(B^)f  and  bending  stiffnesses  (D^)  the  plate  all  depend  upon  the  "governing” 

neutral-surface  position.  Thus,  it  is  possible  to  obtain  closed-form  solutions 
only  for  problems  in  which  the  governing  neutral-surface  positions  remain  constant 


throughout  the  plate.  In  the  case  of  a  rectangular  plar.form,  this  requires  sinu¬ 
soidal  distributions  of  midplane  temperature,  temperature  gradient,  and  pressure. 

It  was  shown  by  Bert  et  al.  (1981a)  that  superposition  is  not  valid  for  blrr.odular 
plates,  so  that  the  usual  Fourler-series  approach  so  popular  in  plate  analysis 
(Timoshenko  and  Woinnwsky-Krieger ,  1959)  is  not  valid.  Fortunately,  however,  the 
finite-element  formulation  is  not  subject  to  this  limitation,  i.e.,  the  neutral- 
surface  positions  may  vary  from  point  to  point  (x,y  position)  as  required.  Thus, 
Reddy  and  Chao  (1980)  considered  rectangular  plates  for  the  more  practical  case  of 
uniformly  distributed  pressure  loading,  and  Reddy  et  al.  (1981)  uniform  distribu¬ 
tions  of  m'.dplane  temperature  and  temperature  gradient.  Reddy  and  Bert  (1932) 
presented  numerical  results  for  clamped  rectangular  plates  subjected  to  uniformly 
distributed  pressure  or  to  point  loading  -it  the  center. 

The  first  dynamic  analyses  of  bimodular  composite-material  plates  were  due  to  Bert 
et  al.  (1981b).  They  presented  both  closed-form  and  finite-element  analyses  of 
rectangular  plates  undergoing  free  vibration.  In  the  case  of  a  single-layer  plate, 
the  plate  exhibits  the  same  stiffnesses  regardless  of  the  sign  of  the  bending 
curvature  in  the  fiber  direction.  Thus,  both  the  positive-curvature  portion  and 
the  negative-curvature  portionofa  cycle  take  place  in  the  same  amount  of  time. 
However,  in  the  case  of  a  two-layer  laminate,  due  to  its  antisymmetric  nature,  the 
plate  stiffnesses  are  different  in  the  two  portions  of  a  cycle  and  the  two  portions 
of  a  cycle  take  place  in  different  lengths  of  time.  Since  the  two  time  intervals 
must  add  up  to  the  total  period  of  one  complete  cycle,  the  frequency  (»i)  defined 
as  2tt/period  is  found  to  be  related  to  the  frequencies  (j|,w2)  of  the  two  portions 
of  a  cycle  by  the  following  simple  relationship : 

u)-1  *  (l/2)(ui-*  +  aj'1)  (ID 

Furthermore,  Bert  eC  al.  (1981b)  demonstrated  that  energy  is  conserved  as  a  plate 
goes  from  the  first  portion  of  a  cycle  to  the  second  portion. 

Recently,  Reddy  (1981)  presented  both  closed-form  and  finite-element  analyses  of 
the  transient  response  of  rectangular  plates  to  a  suddenly  applied,  sinusoidally 
distributed  pressure  loading.  The  linear  matrix  differential  equations  in  time 
were  solved  by  the  well-known  beta  method  due  to  Newmark  (1959).  As  before, 
excellent  agreement  was  achieved  between  the  two  different  types  of  analyses. 

All  of  the  analyses  mentioned  so  far  have  been  linear,  although  piecewise  linear 
(bimodular).  It  is  believed  that  the  first  nonlinear  analysis  of  bimodular 
composite-material  places  is  the  recent  work  of  Turvey  (1981).  He  considered  the 
large  static  deflections  (geometric  nonlinearity)  of  cross-ply,  rectangular,  thin 
plates  using  the  dynamic  relaxation  method  of  solution.  However,  the  numerical 
results  that  he  presented  were  for  material  ehac  was  only  slightly  bimodular  and 
slightly  orthotropic  (e£/e£  =  e£/e£  =  2,  vs.  e£/e£  =  9  to  20  and  E^/E^.  -  100  for 

aramid-rubber,  Bert  and  Kumar,  1931).  Furthermore,  he  did  not  explicitly  state 
which  criterion  for  bimodularitv  (material  model)  he  used.  Recently,  Reddy  and 
Chao  (1982)  independently  considered  the  same  problem,  but  with  thickness  shear 
deformation  included.  They  used  the  finite-element  method  of  solution.  Some 
typical  results  are  shown  In  Fig.  3. 


ORTHOTROPIC  AND  LAMINATED  SHELLS 

The  first  formulation  and  solution  of  a  shell  problem  for  orthotropic  bimodular 
material  was  due  to  Jones  (1971).  He  considered  buckling  of  laminated,  circular 
cylindrical  thin  shells  with  eccentrically  located  stiffeners  and  subjected  to 
arbitrary  combinations  of  axial  and  circumferential  normal  loading.  He  used 
Donnell-type  shell  theory  and  considered  the  stiffeners  to  be  smeared.  Results 


Fig.  3.  Load-deflection  curves  for  thin  square  plates 
bioodular  materials  (a/h  -  100). 

were  presented  in  graphical  form  for  a  ring-stiffened,  two-layer,  c  iv  shell. 

Due  to  the  material  model  used  (Ambartsumyan,  1959),  there  were  cert..,;  discon¬ 
tinuities  on  the  buckling  diagrams.  This  work  was  followed  by  that  of  Jones  and 
Morgan  (1975a),  who  presented  results  for  single-layer  and  two-layer,  cross-ply 
unstiffened  shells.  In  this  work  they  used  two  different  material  models,  the 
Jones  (1977a)  weighted-compliance  model  and  the  Isabekyan  and  Khachatryan  (1969) 
restricted-compliance  model.  Jones  and  Mcigan  (1975b)  also  compared  the  results 
of  their  analyses  to  the  experimental  results  obtained  by  Lenoe-  et  al.  (1969)  on 
three-dimensionally  reinforced  carbon-carbon  cylinders. 

Static  deflection  of  thin  shells  of  orthotropic  bimodular  material  was  analyzed  by 
Karaiya  (1976b).  He  considered  single-layer  circular  cylindrical  shells  subjected 
to  axisymmetric  loading.  Both  cylindrically  curved  panels  and  complete  circular- 
cylindrical  shells  of  single-layer  and  two-layer,  cross-ply  constructed  were 
analyzed  by  Bert  and  Reddy  (1980). 

The  first  treatment  of  orthotropic  bimodular  shells  including  thickness  shear 
deformation  was  due  to  Hsu  et  al.  (1981).  They  treated  circular  cylindrical  shell 
under  both  pressure  loading  and  temperature  variations.  They  introduced  a 
generalized  theory  of  thickness-shear-deformation  shell  theory  (TSDST)  which,  by 
means  of  tracer  coefficients,  reduces  to  the  TSDST  versions  of  the  Sanders  (1959), 
Love  first-approximation  (1927),  and  Donnell  (1933)  thin  shell  theories.  For 
selected  combinations  of  loading  and  boundary  conditions,  both  closed-form  and 
finite-element  solutions  were  presented  for  both  single-layer  and  two-layer,  cross 
ply  constructions.  See  also  Reddy  et  al.  (1981). 

The  first  dynamic  analysis  of  bimodular  shells  was  due  to  Bert  and  Kumar  (1982), 
who  analyzed  the  free  vibration  of  both  thin  and  thick  circular  cylindrical  shells 
of  cylindrical  panel  and  complete  shell  configurations. 


CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 


It  is  concluded  that  many  materials,  especially  fiber-reinforced  composite 
materials,  have  different  stress-strain  behavior  in  tension  and  compression. 
Through  the  bimodular  approximation,  anisotropic  constitutive  relations  are 
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available  and  have  been  successfully  applied  to  a  variety  of  structural  members 
(especially  beams,  plates,  and  shells)  by  means  of  both  classical  and  finite- 
element  methods. 

There  appears  to  be  considerable  need  for  further  research  in  this  area,  mainly: 

1.  Development  of  smooth  nonlinear  (rather  than  bimodular)  con¬ 
stitutive  relations  for  materials  having  drastically  different  behavior  in  tension 
and  compression.  For  analysis  purposes,  these  can  be  treated  as  multimodular, 
i.e.,  piecewise  linear.  Work  on  this  is  currently  in  progress  (Bert  and 
Gordanine j ad ,  19S2). 

2.  Research  on  the  buckling  of  both  columns  and  plates  constructed  of 
bimodular  composite  materials  with  significant  transverse  shear  flexibility. 

3.  Analysis  of  plane  elasticity  problems  such  as  stress  concet'.trat  icn 
at  holes,  fracture  problems,  etc. 

4.  Analysis  of  nonlinear  transient  behavior  of  beams,  plates,  and 
shells  constructed  of  bimodular  composite  materials. 

5.  Analysis  of  edge  effects  at  free  straight  side  edges  and  cutouts 
in  laminates  of  bimodular  composite  materials. 

6.  Analysis  of  infinite  media  of  anisotropic  bimodular  materials  (an 
example  is  the  recent  wave  propagation  analysis  by  Benveniste  (1982). 
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